Quasiparticle relaxation of superconducting qubits in the presence of flux 
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Quasiparticle tunneling across a Josephson junction sets a limit for the lifetime of a superconduct- 
ing qubit state. We develop a general theory of the corresponding decay rate in a qubit controlled 
by a magnetic flux. The flux affects quasiparticles tunneling amplitudes, thus making the decay 
rate flux-dependent. The theory is applicable for an arbitrary quasiparticle distribution. It pro- 
vides estimates for the rates in practically important quantum circuits and also offers a new way of 
measuring the phase-dependent admittance of a Josephson junction. 
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Long coherence times of superconducting qubits rely 
on the decoupling of the order parameter quantum oscil- 
lations from other low-energy degrees of freedom. Quasi- 
particles provide an intrinsic set of states capable of ex- 
changing energy with the qubit degree of freedom. In 
equilibrium, quasiparticle populations should get com- 
pletely depleted at low temperatures, rendering the cor- 
responding qubit relaxation channel ineffective. In prac- 
tice, however, some nonequilibrium quasiparticles are ob- 
served P, Q • The question whether the relaxation driven 
by quasiparticles is comparable to the extrinsic mecha- 
nisms of qubit relaxation has remained open. 

The theory of quasiparticle relaxation was addressed 
in for a charge qubit, whose computational space con- 
sists of two values of charge of the same parity (even or 
odd) stored in a Cooper pair box. The elementary pro- 
cess of relaxation essentially amounts to the well-studied 
quasiparticle poisoning Q,!^: a quasiparticle entering the 
Cooper pair box changes the parity of the state. Later, 
the theory Q was modified to estimate the effect of quasi- 
particles in a "transmon" device @, where the domi- 
nant energy scale comes from the Josephson inductance. 
Quantum fluctuations of phase in a transmon are rela- 
tively small, while the uncertainty of charge in the qubit 
states is significant. The advantage of the transmon is 
its low sensitivity to charge noise. Further consideration 
of the relaxation induced by quasiparticles in supercon- 
ducting qubits was developed in Ref. [2. The properties 
of qubits such as the phase and flux qubits , the trans- 
mon, and the newly developed fluxonium Q can be tuned 
by threading a magnetic flux through the device. 

Here, we predict that the relaxation rates induced 
by quasiparticle tunneling depend on the magnetic flux 
threading the qubit. The theory we develop allows 
for any quasiparticle distribution and any magnitude of 
quantum fluctuations of the phase of the order parame- 
ter; thus, it provides relaxation rates for the entire spec- 
trum of superconducting qubits. We also show that a 
spectroscopic measurement on a device designed to have 
small phase fluctuations may enable one to measure the 
enigmatic phase-dependence of the Josephson junction 



admittance [1]. 

We consider the system consisting of a Josephson junc- 
tion closed by an inductive loop, see Fig.[TJ The Hamilto- 
nian H governing the low-energy dynamics of the system 
can be divided into three parts 
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The first term takes the form of the inductively shunted 
Josephson junction 



= AEc [N - rig] - E,j cos (,5 - 
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where TV = —id/dip is the number operator of Cooper 
pairs passed across the junction, Ug is the dimensionless 
gate voltage, $e is the external flux threading the loop, 
and $0 = h/2e is the flux quantum. With appropri- 
ate choices for the parameters characterizing the qubit 
- charging energy Ec, Josephson energy Ej, and induc- 
tive energy El - and of the biases Ug and $e, Eq- @ 
can describe a single-junction qubit and also a multi- 
junction one, so long as an array of junctions with ener- 
gies Ej Ej is treated as an effective inductance jlol |. 

The quasiparticle term iJqp is the sum of the BCS 
quasiparticle Hamiltonians for the left and right leads 
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FIG. 1: (a) Schematic representation of a qubit controlled by 
a magnetic flux, see Eq. ([2]). (b) Effective circuit diagram 
with three parallel elements - capacitor, Josephson junction, 
and inductor - characterized by their respective admittances. 
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Here a^nai^na) is the quasiparticle annihilation (creation) 
operator, a =t, i accounts for spin, and the quasiparticle 



energies are e->^ = ^ (Cn)^ + {^-'Y^ with and A-' being 
the single-particle energy level n in the normal state of 
lead J, and the gap parameter in that lead, respectively. 
Finally, the tunneling term Ht describes quasiparticle 
tunneling across the junction. 



Ht = t 



H.c. (4) 



The electron tunneling amplitude t here is related to the 
junction conductance g — Ane^v^v^t^ /%. We work in 
the tunneling limit t <^ 1 and assume identical densities 
of states per spin direction in the leads, = = vq. 
The Bogoliubov amplitudes u;^, can be taken real, 
since Eq. (jj]) already accounts explicitly for the phases of 
the order parameters in the leads via the gauge-invariant 
phase difference [§| in the exponentials. Accounting for 
Josephson effect and quasiparticles dynamics by Eqs. ([2])- 
dl]) is possible as long as the qubit energies and charac- 
teristic energies of quasiparticles (measured from A) are 
small compared to A [3|. In this low-energy limit, we 
may approximate ~ i;^ ~ 1 / a/2. Then the operators 
g±Jv/2 ]T;q_ gj) -^^iiich describe transfer of charge ±e 
across the junction, combine to sin((^/2). The superpo- 
sition of the tunneling amplitudes containing e'"^/^ and 
g-«v/2 jg manifestation of interference between tunnel- 
ing of particle- and hole-like excitations, possible in the 
presence of the Cooper pair condensate. We stress that 
the phase (p is an operator subject to quantum fluctua- 
tions, as determined by H^p, not an externally controlled 
parameter as in the "classical" Josephson junction Q. 
Moreover, the non-linear nature of phase-quasiparticle 
coupling is in stark contrast with the linear coupling be- 
tween phase and the electromagnetic environment. 

In general, the qubit Hamiltonian in Eq. ([2]) has a dis- 
crete low-energy spectrum. The tunneling term, Eq. (j4]) , 
couples the qubit to the quasiparticles; therefore, a tran- 
sition between initial, |?), and final, |/), qubit states be- 
comes possible when a quasiparticle is excited during a 
tunneling event. The transition rate Ti^f can be calcu- 
lated using Fermi's Golden Rule, 
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where hujif is the energy difference between the qubit 
states, i?r),qp and i?A,qp are the total energies of the 
quasiparticles in their respective initial {?7}qp and fi- 
nal {A}qp states; double angular brackets ((. . .))qp de- 
note averaging over the initial quasiparticle states whose 
occupation is determined by the distribution functions 

Pm = ((«^^'„t))qp = ((«ni"H))qP = L,R) as- 
sumed to be independent of spin. In terms of the matrix 



elements of the operator sin((^/2) we find 
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In deriving Eq. (jG)) we have made a simplifying assump- 
tion of equal gaps, Al = Ai? = A. More importantly, 
we concentrate on the case of low-lying excitations, as- 
suming the characteristic energy 5E of the quasiparti- 
cles (determined by the distribution functions; in ther- 
mal equilibrium, 5E — ksT) and the energy of the qubit 
transition fiuJif are small, fioJif , 5E <C 2 A. That enables 
us to factorize the transition rate F^.^^ into a product 
of terms which depend separately on the qubit dynam- 
ics and quasiparticle kinetics. The latter determines the 
normalized quasiparticle current spectral density S'qp, 
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X (1 - f§ ((1 + x)A + huj)) +{L^ R) , (7) 

where w > [ll| and we used the relation Ej = gA/Sgx] 
gx = e^/h is the conductance quantum. The inte- 
grand in Sqp equals, up to a factor, the rate of tran- 
sitions of a quasiparticle between the initial and final 
states, properly weighted with their occupation probabil- 
ities (their energies are (1 -I- x)A and (1 -I- a;) A -I- re- 
spectively) . These probabilities are expressed in terms of 
the quasiparticle energy distribution functions, /|;(e) = 

iPiO + /2, with e = ^^eT^ and j - L, R. 

We notice that Sqp is related to the realpart of the 
quasiparticle contribution to the "classical" [9| Josephson 
junction admittance at zero phase difference, 

Sqp(tj) - Sqp(-W) = -— Rerqp(L^), (8) 

71" gK 

at arbitrary ratio fiu/SE and for any quasiparticle dis- 
tribution function. In the rest of this Letter we con- 
sider the "high-frequency" limit 12| Tluj/SE ^ 1 (but 



still huj <C 2A), and equal populations on the two sides 
of the junction, — f^. Then we can simplify Eq. ([5]), 

5qp(w) = -— Reyqp(w) , w > , (9) 

Ti" gK 



(5) and express ReKqp in terms of the quasiparticle density 
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(written using the same approximations as above) as 



Reyqp(a;) = -Xqpg 
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Here Xqp is the quasiparticle density normalized by the 
density of Cooper pairs; in thermal equilibrium, Xqp = 
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yj2'KkBT / IS. f-^l^sT _ Under the above stated assump- 
tions, Eq. pT|) applies to an otherwise arbitrary distribu- 
tion function, and hence to nonequihbrium conditions. 

The structure of S'qp is identical to the corresponding 
one in the Mattis-Bardeen formula [3l , when generalized 
to the case of non-equilibrium distributions [iJl. The 
connection between the Mattis-Bardeen theory, which 
describes absorption of electromagnetic waves imping- 
ing on the surface of a "dirty" superconductor, and the 
theory of qubit relaxation caused by quasiparticles in a 
Josephson junction was noticed in [2]. At the micro- 
scopic level, dissipation in the two settings is caused 
by quasiparticle transitions occurring without momen- 
tum conservation (both systems lack translational invari- 
ance). The difference between the two problems is in 
the form of the perturbation causing the transitions. In 
our case, the quasiparticle transitions are due to the cou- 
pling to the qubit degree of freedom. The matrix element 
(/I sin(i^/2)|i) in Eq. ([B]) characterizes that coupling and 
it is sensitive to flux [15]. 

We stress that Eq. ([6]) holds for any single-junction 
qubit, its properties being encoded in the wavefunctions 
I/) entering the matrix element. For a qubit compris- 
ing multiple junctions the tunneling Hamiltonian Eq. (|4]) 
and hence Eq. ([6]) are given by a sum over junctions, 
with (p being replaced by the phase difference across each 
junction. In a multi-junction qubit, the states \i) and |/) 
depend on all the phase differences (up to the constraint 
set by fluxoid quantization [161]). 

We focus first on the case of a weakly anharmonic sys- 
tem, which already reveals a non-trivial flux dependence 
of relaxation. Its low-lying states, as the example of the 
transmon shows @, can be used as qubit states. We 
assume -El 7^ to eliminate Ug by a gauge transforma- 



tion jl7| . In the transmon Ei^ = 0, but the rig-dependent 



corrections are exponentially small [6|; hence, the results 
below can be applied to a single-junction transmon by 
setting (^0 = 0. In the limit Ec < Ej,El, Eq. ^ de- 
scribes an ideal LC circuit with a junction in parallel. 
The phase across the junction is determined, up to small 
fluctuations, by the external magnetic flux: by minimiz- 
ing the potential energy part of Eq. ^ , we find that the 
phase ipQ at the minimum satisfies 



Ejsin(p^ + EL(ipQ - 27r$e/$o) ^ 0. 



(12) 



Near this minimum, the system behaves as a weakly an- 
harmonic oscillator. Anharmonicity and quality factor Q 
determine the operability of the system as a qubit ^ . In 
the present case, the operability condition can be writ- 
ten as Q/uyj ^ 1, where is the number of levels in 
the anharmonic well. For large Q the system can there- 
fore be operated as a qubit despite its weak anharmonic- 
ity. To the leading order in Ec/{El + Ej cos ipo) [l3] 
and for low- lying levels n <^ n^,, we can neglect anhar- 
monic corrections when evaluating the matrix element in 



Eq. ([6|) [13 . Using standard results for the harmonic os- 
cillator, its value for transitions between two neighboring 
levels is approximated as: 



(n — 1| sin '^W) 
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with Ldio = y^8Ec{EL + Ej cosipo)/h. In deriving 
Eq. (|T^. we replaced sin((^/2) by its linear expan- 
sion, which sets the limit on the excitation level, n <^ 
hujio/Ec- For higher n quantum fluctuations of phase 
are significant and transitions between distant levels pro- 
liferate, due to the nonlinearity of the operator sin(i^/2). 

Concentrating on n ^ Tiujiq/Ec, we use Eq = e^/2C 
and Eqs. ©, dH), and (US]) to find 



n 1 
-1 = — Reyqp(wio) — 
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For the transition between the two lowest states (n = 
1 — > 0) and in the absence of magnetic flux {(fQ = 0), 
Eq. (fT4|) reduces to the expression for the rate of Ref. l2|. 

Using now Eq. (|14|) with n — \, we find the inverse of 
the transition Q-factor at arbitrary </?o, 
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Note that (/Jq and hence the transition frequency ujiq de- 
pend on the external flux <I>e, see Eq. p2)) . We stress 
that the flux dependence presented here is specific to the 
single-junction case; results for multiple-junction qubits 
obtained with the developed method will be presented 
elsewhere. In the limit Ej ^ El, the flux dependence 
can be neglected and the last factor on the right hand 
side of Eq. (fT5)) reduces to « 1. Qubit Q-factors mea- 
sured at temperatures ~ 20 mK are in the range 10^- 
10^ [13] ■ Using typical parameters for Al-based qubits 
(A ~ 2 X 10"" eV, a;io/27r ~ 1 - 10 GHz) we flnd that 
to reproduce the experimental Q-factors with Eq. (fT5)) 
and a thermal equilibrium-like quasiparticle density, we 
must assume an effective quasiparticle temperature an 
order of magnitude larger than the base temperature. 
This points either to the quasiparticles being in a non- 
equilibrium state whose origin is at present unclear, or 
to the prevalence of extrinsic relaxation mechanisms. 

Beside causing dissipation, quasiparticle tunneling 
leads to a shift in the resonant frequency of the cir- 
cuit. In the regime under consideration, the resonant 
frequency is the zero of the total admittance Y{ijj) which 
for the parallel elements representing the qubit (Fig. |T|d) 
isY = Yc + Yl+Yj, with Yc = iwC, Yl = l/iuL, and 
the junction admittance Yj being the sum [§] of a purely 
inductive Josephson term and a quasiparticle term. 



(1 - 2x^^) 

Yj{uj,ip) = — ii-cos(p-f Iqp(w) 
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Here Lj = h/irg/S. is the conventionally defined Joseph- 
son junction inductance, x^^ = fsi^) has the meaning 
of the population of Andreev levels [ilj (x^p = g-A/fceT 
in thermal equilibrium), and 
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is the quasiparticle admittance at zero phase difference 
in the high-frequency limit [its real part agrees with 
Eq. pT|) ]. Both free (xqp) and bound (x^p) quasiparticles 
affect Yj. The frequency shift they cause, measured from 
the frequency in the absence of quasiparticles, is found 
by solving Yc + Yl + Yj 



at linear order in Xqp, x^^: 
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ouj = ^Yqp{ujio) — 
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The imaginary part of frequency shift reproduces - when 
Eq. pT|) for the admittance is used - half the dissipation 
rate Ti^o calculated above, while the real part in the 
high-frequency limit equals 



Re<5a; = i^ 
2 wio 



Xqp (1 - COS(po) 
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where ujp = \/8EjEc /h is the junction plasma frequency. 
The relation Eq. ()18p for duj in the weakly anharmonic 
regime opens new ways to study experimentally the effect 
of quasiparticles on the Josephson junction admittance. 
It may elucidate the phase dependence of its dissipa- 
tive part, where experimental data are still controversial. 
Measuring ImSoj along with KeSuj, Eqs. ([T8l) - ([T9| . may 
also shed light on the nature of nonequilibrium quasi- 
particle distributions: unlike the dissipative part. Re Suj 
depends on both x^^ and Xqp. Experimental efforts to 
verify the flux and quasiparticle density dependence of 
6uj are in progress in our group and elsewhere 22 1. 

The dependence of relaxation on flux is also sensitive 
to the states involved in the transition. The interaction 
between phase and quasiparticles is non-linear (the form 
of the non-linearity is associated with the discreteness of 
charge, as discussed above), thus allowing transitions be- 
tween distant levels even in a harmonic-oscillator poten- 
tial. At Ec/UoJio ^ 1 these transitions are suppressed 
by the smallness of quantum fluctuations of the phase. 
For example, r2^o appears only in the second order in 
fluctuations of (p around ipo, 

r..o = ^-Rer,p(2.,„)f^V^^T^- (20) 
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Here ReFqp is given by Eq. dTTI) . Notice the difference 
in the cpo dependence between Eqs. and (|20p . 



The effect of nonlinear in ip coupling of the qubit de- 
gree of freedom to quasiparticles is striking for a system 
with Ej > Ej^ biased near half the flux quantum [23| . In 
the case of a flux qubit, the potential has a pronounced 
double- well shape and the qubit states |0) and |1) are the 
lowest tunnel-split eigenstates in this potential The 
rate Ti^q vanishes at $e = $o/2 due to the destruc- 
tive interference: at that point, the potential in Eq. ([2]) 
and function sin(^/2 in Eq. ([6]) are symmetric around 
If — TT, while the qubit states |0), |1) are symmetric and 
antisymmetric, respectively. (The symmetry and its con- 
sequences are missed if the sin (p/2 interaction is replaced 
with the linear phase-quasiparticle coupling accepted in 
the environmental approach.) To evaluate ri_j.o at fi- 
nite $e — $0/2 for a qubit governed by Eq. ^ - i.e., a 
single junction connected to an (effective) inductor - we 
consider the case of tunnel splitting [il, 



(21) 



small compared to inductive and plasma energies, hujp 3> 
El 3> eo- Then we can evaluate the matrix element in 
Eq. (jni) in a tight-binding approximation, with the qubit 
states given by linear combinations of states localized in 
the two wells. Using Eq. ([9]), we arrive at 
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where the transition frequency is related to the flux by 



^10 = V^o + [(27r)2SL($e/$o - l/2)r/h and ReFqp is 
given in Eq. (ITT|) . 

Losses in other elements of the qubit give additional 
contributions to the relaxation rate in Eq. (P^ . For ex- 
ample, if the qubit comprises additional Josephson junc- 
tions, by generalizing Eq. (jB)) as previously described one 
can account for the quasiparticle losses in them. Unlike 
Eq. (|22|) those losses remain finite at $e = 4'o/2, but can 
be made small in a chain of iV ^ 1 junctions of larger Ej, 
as they scale with 1/N due to small phase fluctuations 
in the chain (see also (lo|). 

In summary, we have presented a general approach to 
study the quasiparticle relaxation mechanism of super- 
conducting qubits. It enables us to determine how the 
qubit decay rate depends on the magnetic flux used to 
tune the system properties. The method is applicable 
to any superconducting qubit and arbitrary quasiparti- 
cle population. For small phase fluctuations and transi- 
tions between neighboring levels, the decay rate can be 
expressed in terms of the real part of the "classical" Q 
junction admittance, see Eq. (|14p . while the imaginary 
part of the admittance determines the shift in the qubit 
frequency [Eq. (HH)]. This limit is applicable, e.g., to 
a single-junction transmon. Decay rate and frequency 
shift have distinct dependencies on flux and quasipar- 
ticle distribution, so comparing the two quantities may 



give information on the (possibly) nonequilibrium state 
of quasiparticles. 
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